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"Abstract 
-G ■ 

present a new method for the reconstruction of the longitudinal profile of extensive air showers induced by ultra-high energy 
^ Cosmic rays. In contrast to the typically considered shower size profile, this method employs directly the ionization energy deposit 
j_i of the shower particles in the atmosphere. Due to universality of the energy spectra of electrons and positrons, both fluorescence 
and Cherenkov light can be used simultaneously as signal to infer the shower profile from the detected light. The method is 
C$ based on an analytic least-square solution for the estimation of the shower profile from the observed light signal. Furthermore, the 
~ ~ ^extrapolation of the observed part of the profile with a Gaisser-Hillas function is discussed and the total statistical uncertainty of 
shower parameters like total energy and shower maximum is calculated. 



> 

On 
O 

m 
o 

oo 
O 



X 



Key words: Cosmic rays, extensive air showers, air shower reconstruction, air fluorescence, Cherenkov light 
PACS: 96.50.sd 



1. Introduction 

The particles of an extensive air shower excite nitrogen 
molecules in the atmosphere, which subsequently radiate 
ultraviolet fluorescence light isotropically. This fluorescence 
light signal can be measured with appropriate optical de- 
tectors such as the fluorescence telescopes of HiRes [1] , the 
Pierre Auger Observatory [2] or the Telescope Array [3] . 

The number of emitted fluorescence photons is expected 
to be proportional to the energy deposited by the shower 
particles. Recent measurements of the fluorescence yield in 
the laboratory confirm this expectation within the experi- 
mental uncertainties }4|5|6j . Non-radiative processes of ni- 
trogen molecule de-excitation lead to a temperature, pres- 
sure and humidity dependence of the fluorescence yield (see 
e.g. [?])■ For atmospheric parameters of relevance to the 
reconstruction of air showers of ultra-high energy cosmic 
rays, the pressure dependence of the ionization energy de- 
posit per meter track length of a charged particle is almost 
perfectly canceled by the pressure dependence of the fluo- 
rescence yield (see, for example, [8]). Therefore only a weak 
pressure and temperature dependence has to be taken into 
account if the number of emitted photons is converted to 
a number of charged particle times track length, as has 
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been done in the pioneering Fly's Eye experiment j)J . The 
reconstructed longitudinal shower profile is then given by 
the number of charged particles as function of atmospheric 
depth. 

The approximation of assuming a certain number of flu- 
orescence photons per meter of charged particle track and 
the corresponding expression of the longitudinal shower 
development in terms of shower size are characterized by 
a number of conceptual shortcomings. Firstly the energy 
spectrum of particles in an air shower changes in the course 
of its development. A different rate of fluorescence pho- 
tons per charged particle has to be assumed for early and 
late stages of shower development as the ionization en- 
ergy deposit depends on the particle energy [10] . Secondly 
the tracks of low-energy particles are not parallel to the 
shower axis leading to another correction that has to be ap- 
plied [llj . Thirdly the quantity "shower size" is not suited 
to a precise comparison of measurements with theoretical 
predictions. In air shower simulations, shower size is de- 
fined as the number of charged particles above a given en- 
ergy threshold E cut that cross a plane perpendicular to the 
shower axis. Setting this threshold very low to calculate 
the shower size with an accuracy of ~ 1% leads to very 
large simulation times as the number of photons diverges 
for _E cut — > 0. Moreover, the shower size reconstructed from 
data depends on simulations itself since the shower size is 



not directly related to the fluorescence light signal. 

These conceptual problems can be avoided by directly 
using energy deposit as the primary quantity for shower 
profile reconstruction as well as comparing experimental 
data with theoretical predictions. Due to the proportional- 
ity of the number of fluorescence photons to the energy de- 
posit, shower simulations are not needed to reconstruct the 
total energy deposit at a given depth in the atmosphere. 
Another advantage is that the calorimetric energy of the 
shower is directly given by the integral of the energy de- 
posit profile p~2] ■ Furthermore the energy deposit profile 
is a well-defined quantity that can be calculated straight- 
forwardly in Monte Carlo simulations and does not depend 
on the simulation threshold [13] . 

Most of the charged shower particles travel faster 
than the speed of light in air, leading to the emission of 
Cherenkov light. Thus, in general, the optical signal of an 
air shower consists of both fluorescence and Cherenkov 
light contributions. In the traditional method [9] for the 
reconstruction of the longitudinal shower development, the 
Cherenkov light is iteratively subtracted from the measured 
total light. The drawbacks of this method are the lack of 
convergence for events with a large amount of Cherenkov 
light and the difficulty of propagating the uncertainty of 
the subtracted signal to the reconstructed shower profile. 

An alternative procedure, used in |14j . is to assume a 
functional form for the longitudinal development of the 
shower, calculate the corresponding light emission and vary 
the parameters of the shower curve until a satisfactory 
agreement with the observed light at the detector is ob- 
tained. Whereas in this scheme the convergence problems 
of the aforementioned method are avoided, its major dis- 
advantage is that it can only be used if the showers indeed 
follow the functional form assumed in the minimization. 

It has been noted in [15] that, due to the universality of 
the energy spectra of the secondary electrons and positrons 
within an air shower, there exists a non-iterative solution 
for the reconstruction of a longitudinal shower profile from 
light detected by fluorescence telescopes. 

Here we will present an analytic least-square solution for 
the estimation of the longitudinal energy deposit profile of 
air showers from the observed light signal, in which both 
fluorescence and Cherenkov light contributions are treated 
as signal. We will also discuss the calculation of the statis- 
tical uncertainty of the shower profile, including bin-to-bin 
correlations. Finally we will introduce a constrained fit to 
the detected shower profile for extrapolating it to the re- 
gions outside the field of view of the fluorescence telescope. 
This constrained fit allows us to always use the full set of 
profile function parameters independent of the quality of 
the detected shower profile. 

2. Fluorescence and Cherenkov Light Signals 

The non-scattered, i.e. directly observed fluorescence 
light emitted at a certain slant depth Xi is measured at 



the detector at a time i;. Given the fluorescence yield Y? 
[411611715] at this point of the atmosphere, the number of 
photons produced at the shower in a slant depth interval 
AXi is 

N*{X i ) = Yfw i AX i . (1) 
Here, Wi denotes the energy deposited per unit depth at 
slant depth Xi (cf. Fig. Q]) and is defined as 

where dEa ep /dV is the energy deposit per unit volume and 
(ip, R, z) are cylinder coordinates with the shower axis at 
R = 0. The distance interval Azi along the shower axis 
is given by the slant depth interval AA,. The fluorescence 
yield Y? is the number of photons expected per unit de- 
posited energy for the atmospheric pressure and tempera- 
ture at slant depth Xi. The photons from Eq. (Q} are dis- 
tributed over a sphere with surface 47rr?, where denotes 
the distance of the detector. Due to atmospheric attenu- 
ation only a fraction Ti of them reach the detector aper- 
ture with area A. Given a light detection efficiency of e, the 
measured fluorescence light flux y\ can be written as 

yl = diY? Wi AX h (3) 
where the abbreviation d, = eTi . — ■> is used. For the sake 
of clarity the wavelength dependence of Y , T and e will be 
disregarded in the following, but discussed later. 
The number of Cherenkov photons emitted at the shower 
is proportional to the number of charged particles above 
the Cherenkov threshold energy. Since the electromagnetic 
component dominates the shower development, the emitted 
Cherenkov light, , can be calculated from 

N^{X i )=Y l G N^AX t , (4) 

where Nf denotes the number of electrons and positrons 
above a certain energy cutoff, which is constant over the full 
shower track and not to be confused with the Cherenkov 
emission energy threshold. Details of the Cherenkov 
light production like these thresholds are included in the 
Cherenkov yield factor Yp [15|18|19|20| . 

Although Cherenkov photons are emitted in a narrow 
cone along the particle direction, they cover a consider- 
able angular range with respect to the shower axis, because 
the charged particles are deflected from the primary par- 
ticle direction due to multiple scattering. Given the frac- 
tion fc(Pi) of Cherenkov photons per solid angle emitted 
at an angle Pi with respect to the shower axis [18120) . the 
light flux at the detector aperture originating from direct 
Cherenkov light is 

y? d = d i fc(P l )Y l c AX l Nf. (5) 

Due to the forward peaked nature of Cherenkov light pro- 
duction, an intense Cherenkov light beam builds up along 
the shower as it traverses the atmosphere (cf. Fig. [TJ. If 
a fraction f s {Pi) of the beam is scattered towards the ob- 
server it can contribute significantly to the total light re- 
ceived at the detector. In a simple one-dimensional model 




the number of photons in the beam at depth Xj is just the 
sum of Cherenkov light produced at all previous depths Xj 
attenuated on the way from Xj to Xj, by Tji : 

i 

^beampg = £ j.. yC Ax . ^ (g) 
1=0 

Similar to the direct contributions, the scattered Cherenkov 
light received at the detector is then 
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(7) 



,,; </,./;■: u^7 ; ,v< aa\.y;. 

Finally, the total light received at the detector at the time 
ti is obtained by adding the scattered and direct light con- 
tributions: 

yi = vl + y? d + y? s - (8) 



3. Analytic Shower Profile Reconstruction 

The aim of the profile reconstruction is to estimate the 
energy deposit and/or electron profile from the light flux 
observed at the detector. At first glance this seems to be 
hopeless, since at each depth there are the two unknown 
variables Wi and Nf, and only one measured quantity, 
namely yi. Since the total energy deposit is just the sum 
of the energy loss of electrons, Wi and Nf are related via 
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Nf / f e {E,Xi)w e (E)dE, 
Jo 



(9) 



where f c (E,Xi) denotes the normalized electron energy 
distribution and w c (E) is the energy loss per unit depth of a 
single electron with energy E. As is shown in [15 19 20], the 
electron energy spectrum f c (E,Xi) is universal in shower 
age Si = 3/(1 + 2X max /Xi), i.e. it does not depend on the 
primary mass or energy, but only on the relative distance to 
the shower maximum, X max . Eq. §§§ can thus be simplified 
to 

Wi = Nf a,. (10) 



where cti is the average energy deposit per unit depth per 
electron at shower age Si. Parameterizations of ai can be 
found in [10120] . With this one-to-one relation (Eq. [TO]) 
between the energy deposit and the number of electrons, 
the shower profile is readily calculable from the equations 
given in the last section. For the solution of the problem, 
it is convenient to rewrite the relation between energy de- 
posit and light at the detector in matrix notation: Let y = 
(j/i, ?/2, ■ ■ ■ 7 Vn) T be the n-component vector (histogram) 
of the measured photon flux at the aperture and w = 
(wi, W2, ■ ■ ■ , io„) T the energy deposit vector at the shower 
track. Using the expression 
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(11) 



the elements of the Cherenkov- fluorescence matrix C can be 
found by a comparison with the coefficients in equations ([3|) , 
© and 0: 



0, 



i < j 



Ci +c u , i =3 
i, i>3\ 



where 



and 



cf = d t (Yi + f c (/3i)Yf/ai) AX, 



(12) 



(13) 



4 = di / S (A) Tji Yf/ aj AX,. (14) 

The solution of Eq. (fTTj) can be obtained by inversion, lead- 
ing to the energy deposit estimator w: 



(15) 



Due to the triangular structure of the Cherenkov- 
fluorescence matrix the inverse can be calculated quickly 
even for matrices with large dimension. As the matrix ele- 
ments in Eq. fTS]) are always > 0, C is never singular. 
The statistical uncertainties of w are obtained by error 
propagation: 



v w = c- 1 v y (c T ) 1 



(16) 
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Fig. 2. Example of a simulated 10 19 eVproton shower. 



It is interesting to note that even if the measurements yi 
are uncorrelated, i.e. their covariance matrix V y is diag- 
onal, the calculated energy loss values Wi are not. This is 
because the light observed during time interval i does not 
solely originate from Wi, but also receives a contribution 
from earlier shower parts Wj, j < i, via the 'Cherenkov 
light beam'. 



4. Wavelength Dependence 

Until now it has been assumed that the shower induces 
light emission at a single wavelength A. In reality, the 
fluorescence yield shows distinct emission peaks and the 
number of Cherenkov photons produced is proportional to 
pr . In addition the wavelength dependence of the detector 
efficiency and the light transmission need to be taken into 
account. Assuming that a binned wavelength distribution 
of the yields is available (Y^ = j^ k ^^Yi(X)dX), the 
above considerations still hold when replacing cf and c|,- 
in Eq. (HJ) by 

5 4 d = AXi ]T d lk (Yl + fcifa) Yg/cti) (17) 

k 

and 

c£ = AXj £ d lk fM) T jik Y? k /a S) (18) 
fe 

where 



The detector efficiency e k and transmission coefficients 
and Tjik are evaluated at the wavelength Afc . 

5. Validation with Air Shower Simulations 

In order to test the performance of the reconstruction 
algorithm we will use in the following simulated fluores- 
cence detector data. For this purpose we generated proton 
air showers with an energy of 10 19 eV with the C0NEX [21] 



event generator. The resulting longitudinal charged parti- 
cle and energy deposit profiles were subsequently fed into 
the atmosphere and detector simulation package [22] of the 
Pierre Auger Observatory. The geometry and profile of the 
events in this simulated data sample was then reconstructed 
within the Auger offline software framework [23 . 
Only events satisfying basic quality selection criteria have 
been used in the analysis. In order to assure a good re- 
construction of the shower geometry, the angular length 
of the shower image on the camera was required to be 
larger than nine degrees. Moreover, we only selected events 
with at least one coincident surface detector tank (so-called 
hybrid geometry reconstruction [24]). Furthermore we re- 
jected under-determined measured longitudinal profiles by 
demanding an observed slant depth length of > 300 g cm~ 2 
and a reconstructed shower maximum within the field of 
view of the detector. 

An example of a simulated event is shown in Fig. [2] illus- 
trating that the shape of the light curve at the detector 
can differ considerably from the one of the energy deposit 
profile due to the scattered Cherenkov light detected at 
late stages of the shower development. The reconstructed 
energy deposit curve, however, shows on average a good 
agreement with the generated profile. 
Since longitudinal air shower profiles exhibit similar shapes 
when transformed from slant depth X to shower age s (see 
for instance [IS]), a good test of the profile reconstruction 
performance is to compare the average generated and re- 
constructed energy deposit profiles as a function of s nor- 
malized to the energy deposit at shower maximum. As can 
be seen in Fig. ([3]), the difference between these averages 
is < 1.5% and it can be concluded that the matrix method 
introduced here performs well in reconstructing air shower 
profiles without a prior assumption about their functional 
shape. 




shower age 

Fig. 3. Average generated and reconstructed energy deposit profiles. 

6. Shower Age Dependence 

Due to the age dependence of the electron spectra 
f e (E, Si), the Cherenkov yield factors Yp and the average 
electron energy deposits a.i depend on the depth of shower 
maximum, which is not known before the profile has been 
reconstructed. Fortunately, these dependencies are small: 
In the age range of importance for the shower profile recon- 
struction (s G [0.8, 1.2]) a varies by only a few percent [20] 
and Y c by less than 15% [15] ■ Therefore, a good estimate 
of a and Y c can be obtained by setting s = 1 over the full 
profile or by estimating A max from the position maximum 
of the detected light profile. After the shower profile has 
been calculated with these estimates, A max can be deter- 
mined from the energy deposit profile and the profile can 
be re-calculated with an updated Cherenkov-fluorescence 
matrix. The convergence of this procedure is shown in 
Fig. El After only one iteration the A max (energy) differs 
by less than 0.1 gcm~ 2 (0.1%) from its asymptotic value. 
Note that age dependent effects of the lateral spread of the 
shower on the image seen at the detector [26127] . though 
not discussed in detail here, have also been included in the 
simulation and reconstruction. 

7. Gaisser-Hillas Fit 

A knowledge of the complete profile is required for the 
calculation of the Cherenkov beam and the shower energy. 
If due to the limited field of view of the detector only a part 
of the profile is observed, an appropriate function for the 
extrapolation to unobserved depths is needed. A possible 
choice is the Gaisser-Hillas function [28] 
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Fig. 5. X max and energy difference with respect to the tenth shower 
age iteration. 

which was found to give a good description of measured 
longitudinal profiles [59]. It has four free parameters: A max , 
the depth where the shower reaches its maximum energy 
deposit w max and two shape parameters A and A. 
The best set of Gaisser-Hillas parameters p can be ob- 
tained by minimizing the error weighted squared difference 
between the vector of function values fcH and w, which is 



X^ H = [w-f(p)] T Vw 1 [w-f(p)] 



(21) 



This minimization works well if a large fraction of the 
shower has been observed below and above the shower max- 
imum. If this is not the case, or even worse, if the shower 
maximum is outside the field of view, the problem is under- 
determined, i.e. the experimental information is not suffi- 
cient to reconstruct all four Gaisser-Hillas parameters. This 
complication can be overcome by constraining An and A to 
their average values (An) and (A). The new minimization 
function is then the modified \ 2 
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Vx 



(22) 



where the variances of An and A around their mean values 
are in the denominators. 

In this way, even if Xgh ^ s n °t sensitive to A and A, the 
minimization will still converge. On the other hand, if the 
measurements have small statistical uncertainties and/or 
cover a wide range in depth, the minimization function is 
flexible enough to allow for shape parameters differing from 
their mean values. These mean values can be determined 
from air shower simulations or, preferably, from high qual- 
ity data profiles which can be reconstructed without con- 
straints. 

Eq. ([22]) can be easily extended to incorporate correlations 
between A and A and the energy dependence of their mean 
values. Air shower simulations indicate a small logarithmic 
energy dependence of the latter (< 25% and 5% per decade 
for An and A respectively [3D])- In practice it is sufficient to 
use energy independent values determined at low energies, 
because at high energies the number of measured points 
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Fig. 4. Energy and X max reconstruction accuracy as function of the amount of detected Cherenkov light. 



is large and thus the constraints do not contribute signifi- 
cantly to the overall x 2 ■ 

The accuracy of the reconstructed energy, obtained by 
integrating over the Gaisser-Hillas function (see below), 
and that of the depth of shower maximum, are displayed 
in Fig. |4] as a function of the relative amount of Cherenkov 
light. Note that the good resolutions of « 7% and 20 gem -2 
are of course not a feature of the reconstruction method 
alone, but depend strongly on the detector performance 
and quality selection. The mean values of difference to the 
true shower parameters, however, which are close to zero for 
both fluorescence and Cherenkov light dominated events, 
indicate that both light sources are equally suited to recon- 
struct the longitudinal development of air showers. 
A slight deterioration of the resolutions can be seen for 
events with a very small Cherenkov contribution of < 10%. 
Such showers are either inclined events which developed 
high in the atmosphere where the light scattering probabil- 
ities are low or deep vertical showers, for which most of the 
late part of the shower is below ground level. Both topolo- 
gies result in a somewhat worse resolution: The former cor- 
respond to larger than average distances to the detector 
and the latter to shorter observed profiles. 



8. Error Propagation 

A realistic estimate of the statistical uncertainties of im- 
portant shower parameters is desired for many purposes, 
like data quality selection cuts, or the comparison between 
independent measurements like the surface and fluores- 
cence detector measurements of the Pierre Auger Obser- 
vatory or the Telescope Array. The uncertainties of u> ma x, 
A max , Ao and A obtained after the minimization of Eq. (|22p , 
reflect only the statistical uncertainty of the light flux, 
which is why these errors will be referred to as 'flux uncer- 
tainties' (uflux) in the following. Additional uncertainties 
arise from the uncertainties on the shower geometry (tT geo ), 
atmosphere (cr a tm) and the correction for invisible energy 

(Cinv)- 



8.1. Flux uncertainty of the calorimetric energy 

Even with the flux uncertainties of the Gaisser-Hillas 
parameters it is not straightforward to calculate the flux 
uncertainty of the calorimetric energy, which is given by 
the integral over the energy deposit profile: 



£ ca i= / Zgh(A) dA 
Jo 



(23) 



To solve this integral one can substitute 

t= A_Ao ^ ^ X max -X (M) 
A A 
in the Gaisser-Hillas function Eq. (|20|) to get 



/GH(t) = W max ( 7 ) e 



(25) 



which can be identified with a Gamma distribution. There- 
fore, the above integral is given by 

£ C ai = Auw T(£ + l), (26) 

where T denotes the Gamma-function. Thus, instead of do- 
ing a tedious error propagation to determine the statistical 
uncertainty of E ca \ one can simply use it directly as a free 
parameter in the fit instead of the conventional factor w max : 



Ar« + i) e ~' 15 < 27 > 

In this way, crfiux(-Ecai) is obtained directly from the x 2 + 1 
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contour of Eq. ([22 



8.2. Geometric uncertainties 

Due to the uncertainties on the shower geometry, the 
distances fj to each shower point are only known within 
a limited precision and correspondingly the energy de- 
posit profile points are uncertain due to the transmission 
factors T(ri) and geometry factors l/(47rr 2 ). Further- 
more, the uncertainty of the shower direction, especially 



the zenith angle 9, affects the slant depth calculation via 
Aslant = X vcrt / cos 9 and thus X max . Finally, the amount 
of direct and scattered Cherenkov light depends on the 
shower geometry, too, via the angles /Jj. 

The algorithms used to reconstruct the shower geome- 
try from fluorescence detector data usually determine the 
following five parameters [9j, irrespective of whether the 
detectors operate in monocular, stereo or hybrid mode: 



a = {0sdp,$sdp,?o, Rp,Xo}- 



(28) 



9s dp and $5 r> p are the angles of the normal vector of a 
plane spanned by the shower axis and the detector (the so 
called shower-detector-plane), xo denotes the angle of the 
shower within this plane and To and R p are the time and 
distance of the shower at its point of closest approach to 
the detector. 

For any function q(a) standard error propagation yields the 
geometric uncertainty 
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(29) 



where V a denotes the covariance matrix of the axis param- 
eters. As the calorimetric energy and X max depend non- 
trivially on the shower geometry, the above derivatives need 
to be calculated numerically, i.e. by repeating the profile 
reconstruction and Gaisser-Hillas fitting for the ten new 
geometries given by a« ± WV? = a, ± <7j to obtain 



Aj = - — Oi 

da,; 



1 



q(ai + ai) - - <Ti)\ 



with which Eq. 



reads as 

5 5 

i=i j=i 



where 



Pij 



V a 



(30) 



(31) 



(32) 



V " 33 

denote the correlation coefficients of the geometry param- 
eters a, and aj . 




Mean 
RMS 



(a) Energy. 



(b) X n 



Fig. 6. Pull distributions 
8.4. Invisible energy 

Not all of the energy of a primary cosmic ray particle ends 
up in the electromagnetic part of an air shower. Neutrinos 
escape undetected and muons need long path lengths to 
fully release their energy. This is usually accounted for by 
multiplying the calorimetric energy, Eq. (|23[) . with a cor- 
rection factor /; nv determined from shower simulations to 
obtain the total primary energy 



finv E CSl \. 



(33) 



The meson decay probabilities, and thus the amount of neu- 
trino and muon production, decrease with energy, therefore 
/i nv depends on the primary energy. For instance, in [33j it 
is parameterized as 

-1 



/ inv = (a + bE c cal Y 



(34) 



where a, b and c denote constants depending on the pri- 
mary composition and interaction model assume cCD. This 
energy dependence needs to be taken into account when 
propagating the calorimetric energy uncertainty to the to- 
tal energy uncertainty. 

Due to the stochastic nature of air showers, the correc- 
tion factor is subject to shower-to-shower fluctuations. The 
statistical uncertainty of /; nv was determined in [34] and 
can be parameterized as follows: 

tr(/inv) « 1-663 • 10 6 • lg(£W eV)" 6 36 . (35) 

Typical values are 2.5% at 10 17 eV and 0.9% at 10 20 eV. 



8.3. Atmospheric uncertainties 



8.5. Total statistical uncertainty 



Whereas the Rayleigh attenuation is a theoretically well 
understood process, the molecular density profiles and 
aerosol content of the atmosphere vary due to environ- 
mental influences and need to be well monitored in order 
to determine the slant depth and transmission coefficients 
needed for the profile reconstruction. Uncertainties in 
these measured atmospheric properties (see for instance 
|31l32j ) can be propagated in the same way as the geo- 
metric uncertainties by determining the one sigma shower 
parameter deviations via Eq. ((30|) . 



Summarizing the above considerations, the statistical 
variance of the total energy is 



cr stat(^tot) - #tot 0- 2 (/inv) 



+ /inv 5> 2 (£ cal ), (36) 



1 Note that here only the statistical uncertainties of the invisible 
energy correction are discussed. For an estimate on the related sys- 
tematic uncertainties see 1341 



v 



where i runs over the geometric, atmospheric and flux un- 
certainties. Since the invisible energy correction does not 
affect the depth of shower maximum, its uncertainty is sim- 
ply given by 



Cstat (^max) — , ] of (^mai) 



(37) 



Again we use simulated events to verify the validity of the 
above considerations. The pull distributions of the recon- 
structed energy and shower maximum, shown in Fig. I37( 
both have a width of approximately one, which means that 
the total uncertainties from Eqs. (f3"6")) and (f3"7| are good 
estimators for the actual event-by-event measurement un- 
certainties. 

9. Conclusions and Outlook 

In this paper a new method for the reconstruction of 
longitudinal air shower profiles was presented. With the 
help of simulations we have shown that the least square 
solution yields robust and unbiased results and that un- 
certainties of shower parameters can be reliably calculated 
for each event. 

Events with a large Chcrcnkov light contribution are cur- 
rently usually rejected during the data analysis (see for 
instance |14I35| ) However, as we have shown, there is no 
justification for rejecting such showers, once experimen- 
tal systematic uncertainties are well understood. Because 
events with a large Cherenkov contribution have different 
systematic uncertainties to those dominated by fluores- 
cence light, both event classes can be compared to study 
their compatibility. 

At energies below 10 17 5 eV, where new projects [36 37 
are planned to study the transition from galactic to ex- 
tragalactic cosmic rays, events with a large fraction of 
direct Cherenkov light will dominate the data samples, 
because the amount of light, and thus trigger probability, 
of these events is much larger than that of a fluorescence 
dominated shower. If at these energies it is still possible to 
measure an accurate shower geometry, the fluorescence de- 
tectors should in fact be used as Cherenkov-Fluorescence 
telescopes. 
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